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Critical behavior of the contact process is studied in annealed scale-free networks by mapping 
it on the random walk problem. We obtain the analytic results for the critical scaling, using 
the event-driven dynamics approach. These results are confirmed by numerical simulations. The 
disorder fluctuation induced by the sampling disorder in annealed networks is also explored. Finally, 
we discuss over a possible discrepancy of the finite-size-scaling theory in annealed and quenched 
networks in spirit of the droplet size scale and the linking disorder fluctuation. 
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I. INTRODUCTION 

Critical phenomena in complex networks have been at- 
tracting a lot of interest [l| . Complex networks are char- 
acterized by a so-called small- world property The 
number of neighbors of a node increases exponentially 
with the distance from it. For this property, it is be- 
lieved that critical phenomena in complex networks be- 
long to the mean field universality class. Nevertheless, 
structural heterogeneity leads to rich behaviors. For ex- 
ample, in scale-free (SF) networks having a power-law 
degree distribution P{k) ^ 3], mean-field critical 

exponents may vary with the degree exponent 7 [l| . 

Recent studies raise an important issue on the finite- 
size-scaling (FSS) theory in complex networks A 
scale-free network with N nodes has a maximum cutoff 
fcmax in degree. In most cases without any constraint, the 
cutoff scales as fc,nax A^i/^T"!), which is determined by 
the condition J2k>k ^(^) = l/N. This is called the 
natural cutoff. One may impose a forced cutoff 



k — AT^/'" 

"-max — J-^ 



(1) 



with the cutoff exponent u > 7 — 1. Taking the ther- 
modynamic limit, one should take the limit N 00 and 
fcmax 00 simultaneously. This may give rise to an in- 
tricate finite-size effect M ■ 

Hong, Ha, and Park [j| developed a FSS theory based 
on the single-parameter scaling hypothesis. Their the- 
ory predicts the values of the FSS exponents in the Ising 
model (including more general equilibrium 0" theory) 
and the contact process (CP), respectively. The CP is a 
reaction-diffusion model describing an epidemic spread- 
ing, which exhibits a prototype nonequilibrium phase 
transition from an inactive phase into an active phase Q . 
It has been suggested that the FSS exponents depend 
only on the degree exponent 7, regardless of the cutoff if it 
is not strong enough (co < 7). Note that this condition in- 
cludes networks with the natural cutoff {lo = 7— 1) as well 
as networks with a weak forced cutoff (7 — 1 < < 7) . 
These results were confirmed numerically in the static 



model 0] having the natural cutoff and the uncorrelated 
configuration model (UCM) [|. 

Castellano and Pastor-Satorras [H] considered the CP 
in the so-called random neighbor (annealed) network. 
Links are not fixed but fluctuate in this annealed net- 
work. At each time step, neighbors of a node are cho- 
sen independently and randomly according to the degree 
distribution. It contrasts with a network where links are 
fixed permanently in time once they are formed. In or- 
der to stress the distinction, the former network will be 
referred to as an annealed network, while the latter net- 
work as a quenched network. From the analysis of the 
survival probability at the critical point, they found that 
the dynamic exponent characterizing the relaxation time 
depends not only on 7 but also on to when w > 7 — 1 
(all networks with a forced cutoff) and 2 < 7 < 3. In 
particular, it has been shown that there are two different 
characteristic time (and also the order parameter) scales 
which make a single-parameter scaling impossible. From 
the relaxation time scaling, the order parameter in the 
quasi-steady state also scales with N with an exponent 
depending on both 7 and lo. 

At a glance, the results of Refs. Q and seem incom- 
patible (single-parameter versus two-parameter scaling 
and cutoff-independent versus cutoff-dependent scaling) 
when 7 — 1 < a; < 7 (weak forced cutoff) and 2 < 7 < 3 
(highly heterogeneous regime). But it is not true. The 
FSS theory of Ref. [J concerns a quenched scale- free net- 
work, while that of Ref. [H] concerns an annealed scale- 
free network. Quenched disorder in linking topology gen- 
erates local correlations through quenched links between 
nodes, which are responsible for the shift of the phase 
transition point and its disorder fluctuations. Therefore, 
one may not rule out a possibility that the disorder fluc- 
tuations near the phase transition point may wipe away 
or at least significantly alter the cutoff-dependent scaling 
regime, see Sec. VI. 

In this paper, we present a full FSS theory govern- 
ing the critical and off-critical scaling behaviors of the 
CP in annealed networks. In Sec. II, an annealed net- 
work is introduced without any sampling disorder and 
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the heterogeneous mean field theory is briefly reviewed 
for the CP. The critical dynamics is analyzed in Sec. Ill, 
while the off-critical scaling is investigated in Sec. IV. In 
Sec. V, we discuss the effect of sampling disorder in an- 
nealed networks and its self-averaging property. Finally, 
we summarize our results along with a brief discussion 
on the effect of linking disorder in quenched networks. 



II. CP IN ANNEALED NETWORKS 

We consider the annealed scale-free networks with the 
degree distribution P{k) = ak~^ for fcmin l£ k < /cmax 
with a normalization constant a and P{k) = elsewhere. 
The maximum degree fcmax scales with network size N 
as in Eq. ((!]) and the minimum degree fcmin is an 0(1) 
constant. Since neighbors of each node need not be speci- 
fied, an annealed network is realized by choosing a degree 
sequence {ki, • • • , k^} only. 

There are two different ways in choosing the de- 
gree sequence. One may assign degree k to Nk nodes 
deterministically in such a way that '^ki^f^Nk' = 
mt[N J2k'>k ^(^')] foi' ^ the decreasing order start- 
ing from /cmax, where int[a;] is the integer part of x. One 
may easily show that the maximum degree realized us- 
ing this assignment algorithm is the same order in A'^ of a 
given fcmax when uj > 7 — 1. Or, one may draw probabilis- 
tically N values of k in accordance with the probability 
distribution P{k). The probabilistic method yields an 
ensemble of different samples, which makes an ensem- 
ble average necessary. We mainly consider the annealed 
network realized by the deterministic method. Sample- 
to-sample fluctuations in the ensemble generated by the 
probabilistic method will be discussed in Sec. V. 

The CP on the annealed SF network is defined as fol- 
lows. Each node is either occupied by a particle or empty. 
A particle on a node is annihilated with probability p 
or branches one offspring to its neighbor, if empty, with 
probability (1 — p). At each time step, a neighbor of a 
node is selected among all other nodes with probabili- 
ties proportional to their degree. Since a node is coupled 
only probabilistically with all other nodes, the mean field 
theory becomes exact in the annealed network. 

Let n{t) be the number of particles at time t. Follow- 
ing Ref. [5[ in a quasistatic approximation for large t, it 
increases by 1 with probability 



w+ 



kP{k) 



(k) l + Xpk/{k)'' 



or decreases by 1 with probability 

w_ = pp 



(2) 



(3) 



after a time step At ~ 1/N. Here p — n/N is the particle 
density and A = {l—p)/p with the mean degree (fc). 

The transition probability w+ contains a nontrivial p- 
dependence. When the thermodynamic limit is taken 



first [3| or the density is high (p » l/fcmax) in finite 
networks one may arrive at a singular expansion 



w+/p = 

with a constant c and 



Xp~cp^-^ + 



i{7,3}. 



(4) 



(5) 



When the density is low (p ^ l/^max) in finite networks, 
one can expand the denominator in Eq. Q to obtain 



w+/p ~ Xp - X^gp^ + 



(6) 



where g = {k^)/{k)^ with {k") = Y,k k"P(k). Note that 
g is an 0(1) constant for 7 > 3, while it scales as 17 ~ 
^max ~ 7V(3-7)/c^ for 2 < 7 < 3. The scaling behavior 
can be rewritten as 



(7) 



for general 7(7^ 3) and ui > j ~ 1. At j — 3, g ^ log N. 

As the stochastic fluctuation (Ap)/p (multiplicative 
diffusive noise) becomes negligible in the N 00 limit, 
one can write the rate equation for the average particle 
density in the continuum limit as 



dp 

dt ^ 



(8) 



It is clear that the system undergoes an absorbing phase 
transition at p = Pc — 1/2 (Ac = 1) at all values of 7 > 2 
in the thermodynamic limit. The particle density near 
the critical point scales as p ~ (A — Ar)^ with the order 
parameter exponent /9 = 1/(0 — 2) At 7 = 3, an 
additional logarithmic correction appears as p ~ 
Ac)/|log(A-Ac)|. 



III. CRITICAL DYNAMICS 

We consider the CP at the critical point {p ~ 1/2 or 
A = 1). One may regard the particle number n (0 < 
n < N) as a coordinate of an one-dimensional random 
walker [5i]. At each time step At = 1/N, the walker 
jumps to the right with probability i/7-|_ or to the left 
with probability W-, or does not move with probability 
1 — (w-i- -|- w_ ) . The walker is bounded by an absorbing 
wall at n = and a reflecting wall at n = A''. Reaching 
the absorbing wall, it will be trapped there forever. 

It turns out that an event-driven dynamics is useful. In 
this dynamics, the walker always makes a jump at each 
time step At = 1 to the right or left with probabilities 



'W± 



'W± 



(9) 



This is equivalent to the original problem if one rescales 
the time with the relation 



dt : 



1 



dr 



N w+ + w- 



(10) 
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A. Defect dynamics 

It is interesting to study how particles spread starting 
from a localized seed. Dynamics initiated from a single 
particle is called the defect dynamics [1, Q . So its initial 
condition is n(0) = uq = 1. 

Quantities of interest are the survival probability Ps (t), 
the probability that the system is still active at time t, 
and Usit), the number of particles averaged over surviv- 
ing samples. At the critical point, they exhibit power-law 
scalings 

Psit) - t-^ and n,(t) ~ (11) 

for t < tc{N) with the relaxation time scaling as 

tc ~ (12) 

At t — tc, the system starts to feel its finite size and 
ns{t) saturates. For t > tc, Ps{t) decays exponentially. 
The critical exponents S, fj, and z are universal. Note 
that fj = 6 + rj where rj is the particle number growing 
exponent for all samples. 

Initially, po = ^q/N is so small (much smaller than 
1/fcmax) that one can always use the expansion in Eq. ^ 
for 'W+ . We will confirm that this is valid throughout the 
defect dynamics. Then, in the event-driven dynamics, 
the jumping probability for the walker at site n is given 

by 



1 



gp 



and 



1 



for small p. This shows that the walker performs biased 
walks toward the absorbing wall with the drift velocity 



Vdrift 



dn 



gp 



(13) 



The bias is negligible {vdrift/w± <ti 1) during the ini- 
tial stage since gp <^ 1. Hence, for sufficiently small r, it 
suffices to consider the unbiased random walk motion in 
the presence of the absorbing wall at n = 0. The effect 
of the absorbing wall can be taken into account by using 
the image method [10]. This yields that the surviving 
probability decays as 



P,(t) ~no(W2) 



-1/2 



(14) 



and that the surviving walker spreads out diffusively as 



,(t) ~ VW2- 



(15) 



by 



The diffusion velocity for the surviving walkers is given 



Vd 



If fuse = ns{T)/T ~ ^7r/(2T) ~ 7r/(2ns). 



(16) 



As r increases, the diffusion velocity becomes smaller 
while the bias becomes stronger. The walker reaches a 
stationary state when the diffusion velocity and the drift 



velocity are balanced. The condition Vdiffuse ~ Ivdri/tl 
yields that the walker reaches the stationary state at po- 
sition 



and at time 



^/Wg 



N/g. 



(17) 



(18) 



This result is self-consistent with the underlying assump- 
tions that pfcmax ^ 1 and gp 1 . 

The time scales t and r are related through Eq. pop . 
Using Eqs. p5|) and (fTO|l . one finds that 



t 



n{T') 

Therefore we conclude that 

Psit)^not-\ 



which leads to 



g, 



S = 1, 

z = (1- (3-6')/tj)/2. 



(19) 



(20) 
(21) 
(22) 

(23) 
(24) 
(25) 



The result for 8 and z coincides with that of Ref. 0. At 
^ = J,^tc^nf ^ (TV/ log iV) 1/2. 



B. Static dynamics 

The static dynamics starts with the initial condition 
that all nodes are occupied, uq — N {po — 1). We con- 
sider the scaling behavior of ps, the particle density av- 
eraged over surviving samples. 

The rate equation ([8]) takes a different form depending 
on the particle density. When pky^ax ^ 1, it becomes 
dp/dt — — cp^~i/2, which yields 



Ps{t)^t-'/(o-^\ 



(26) 



If the density becomes sufficiently small such that 
pfcmax ^ 1, then the rate equation should be replaced 
by dp/dt = —gp^/2, which yields the solution 



Psit) ^ {gt)-\ 



(27) 



The crossover between the two regimes takes place at 
time 



t, ^ g(«-2)/(3-e) ^ AT--, with z, = {e- 2)/t 



(28) 



At this crossover time scale t = t^,, the system starts to 
feel the finite upper bound of the maximum degree, fcmax- 
The density at the crossover is given by 



-1/(3-9) ^ ^-a. ^-^^ 



l/uj. (29) 
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FIG. 1: (a) Schematic plot of ps vs. t in the log-log scale at 
the critical point in the annealed networks with 2 < 7 < 3 
and u) > 7 — 1. The solid (dashed) line corresponds to the 
static (defect) dynamics, (b) Schematic plot of ps vs. e in the 
same condition. 



Finally, the system reaches the stationary state. From 
Eq. (flT)) . the particle density at the stationary state is 
given by 

pf - VTJigN) - iV~" with a ^ (l + (3-6')/w)/2. (30) 

The saturation time tc determined from (gic)^^ ~ iV^" 
has the same scaling behavior as the relaxation time in 
the defect dynamics (see in Eqs. ([22|)and (f25|) ). This 
means that the finite systems reach the stationary state 
at the same time scale, irrespective of the initial condi- 
tions. 

There are a few remarks. For a; > 7 — 1 and 7 < 
3 {Q = 7), there exist two distinct A^-dependent time 
scales ~ 7V^* and tc ^ with < z. The former 
comes into play due to the finiteness of the maximum 
degree k^ax_ ~ N^^'^, while the latter is the time scale 
to reach the stationary state in finite networks. This 
implies that finite-size effects in the annealed SF net- 
works depends on the limiting procedure how N and fcmax 
are taken to infinity. For 7 > 3 (6* = 3), the distinc- 
tion between the first regime (Eq. ([^5]) ) and the second 
regime (Eq. ([27)) ) disappears. The particle density de- 
cays as Ps ~ for t < with z = 1/2, and then 
saturates to the stationary state value pf ^ N~'^/'^. At 
7 = 3, - {tlogt)-^ for t < tc - (iV/logiV)i/2 and 
pf - [NlogN]-^''^. 

The systems with the natural cutoff (a; ~ 7 — 1) are 
special. Even for 7 < 3, the two time scales and 
coincide having z.^, — z — (7 — 2)/ (7 — 1). This means 
that the second regime does not exist. The density decays 
as Ps ^ t~'^/h~^) for t < N^, and then saturates to 
pf - Ar-i/(7-i). 

The defect and static dynamics at criticality are illus- 
trated schematically in Fig.[IJa). 



C. Numerical simulations 

We have performed numerical simulations in the an- 
nealed SF networks to confirm the analytic results. In the 
defect simulations, the survival probability is expected to 
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FIG. 2: (Color online) (a) Scaling plot for the survival prob- 
ability, (b) Scaling plot for the number of particles averaged 
over surviving samples. The upper, middle, and lower sets of 
data correspond to the annealed SF networks with 7 = 2.5 
and ui — 1.5, 2.0, and 3.0, respectively. For readability, each 
data set is scaled down by a constant factor. The dashed lines 
are guides to the eyes having slope —1 in (a) and 1 in (b). 
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FIG. 3: (Color online) (a) Critical density decay at 7 = 2.5 
and Lo = 2.5. The network sizes are TV = 10^, • • • , lO'^. Ef- 
fective exponent q versus t/N^' at 7 = 2.5 and (b) uj — 2.0, 
(c) 2.5, (d) 3.0. 



scale as 



Psit,N) = N-^^Vit/N^) 



(31) 



The scaling function behaves as 7-" (a:) —f and 

decays exponentially as x —^ 00. The particle number 
averaged over surviving samples is expected to scale as 



(32) 
and 



The scaling function behaves as Af{x) ^ x"^' a.s x 
converges to a constant as a; — *■ cxd. 

Figure [5] shows the scaling plots according to the scal- 
ing form in Eqs. pip and p2p with the exponent values in 
Eqs. (1121), (El), and The annealed networks of size 

N — 10^, • • • , 10^ were generated with the deterministic 
method. The plotted data were obtained by averaging 
over 10^ runs. The nice data collapse confirms the valid- 
ity of the analytic result. 
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We have also performed the static simulations. Numer- 
ical data at the critical point obtained in networks with 
7 — 2.5 and uj = 2.5 are presented in Fig. [3ja). Unlike in 
the schematic plot in Fig. [2 the crossover between two 
regimes with ps ~ ^-1/(^-2) ^^^^ ^ ^-i^ respectively, 
is not prominent. Moreover, the decay exponents seem 
to deviate from the expected values significantly. In or- 
der to understand the origin of the discrepancy, we have 
performed a local slope analysis. As an estimate for the 
density decay exponent, we define an elTective exponent 
q{t) = — hi{ps{t) / ps{t /m)) / Inm with a constant m = 4. 
The effective exponents measured a.t u; — 2.0, 2.5, and 
3.0 are plotted in Figs. [3jb), (c), and (d), respectively, 
against the scaling variable x — t/N^* . The analytic the- 
ory predicts that q should converge to 1/(7 — 2) = 2 (1) 
as N increases for x < 1 (x > 1). The effective exponent 
plot shows a weak but clear tendency toward the analytic 
prediction. For x < 1, the effective exponents steadily in- 
creases above 1 with network size, but still much lower 
than the predicted value 2 even for N ^ 10^. Moreover 
there is no appreciable power-law region (flat region for 
q). For a; > 1, the effective exponents overshoot the pre- 
dicted value 1 till N = 10^, but start to decrease slightly 
at iV = 10^. We also notice the appreciable flat region in 
this case. 

In order to identify numerically the power-law scaling 
in the first regime, it would be required that at least 

~ 10^ (two log decades). With 7 = 2.5 and w = 2.5, 
the system size must be larger than ~ 10"'^°, which is 
beyond the current computer capacity. 

We have also studied the FSS behavior of the 
stationary-state particle density at criticality. It exhibits 
a power-law scaling with N as pf^ ^ 7V~". It is found 
(not shown here) that the numerical result for the ex- 
ponent a is compatible with the analytic result given in 
Eq. ([50)1 . However, a discrepancy becomes noticeable as 
7 becomes smaller and co becomes larger, due to strong 
finite size effects. The degree distribution becomes singu- 
lar as 7 approaches 2. At large to, the maximum degree 
fcmax ~ N^/^ grows SO slowly that it becomes difficult to 
observe the asymptotic scaling behavior. 



IV. OFF-CRITICAL SCALING 

For 2 < 7 < 3, the particle density exhibits distinct dy- 
namic characteristics depending on whether yofcmax > 1 or 
pfcmax < 1- This causes an interesting cutoff-dependent 
FSS behavior at the critical point. Such a cutoff de- 
pendence disappears far from the critical point. How- 
ever, in finite systems near the critical point, the cutoff- 
dependence can still survive to lead to an anomalous FSS 
behavior. For 7 > 3, the system shows a simple normal 
FSS behavior. 

Near the critical point at p — Pc(l ^ £), the rate equa- 



tion for the density, Eq. ([8]), reads for 2 < 7 < 3 

dps/dt = eps - d pY^ for Ps^max > 1, (33) 
= eps - \gpi for psk^a^ < 1, (34) 

where fc^ax = N^^'^, g = {k^)/{kf - ^^(3-7)/^^^ and 
c' = c/2. 

By setting dps/dt = 0, one obtains that the stationary 
state solution is given by 

pf ^ for p,A:„iax > 1, (35) 

~ 2e/g for p^Aimax < 1 (36) 

with the bulk order parameter exponent 

/?-l/(7-2). (37) 

This shows that the stationary state solution also de- 
pends on the degree cutoff fcmax- There is a crossover at 
e = e* with 

~ 5-('r-2)/(3-7) _ iY-iZ-^-with = (38) 

For e < the order parameter scaling changes into the 
7-independent ordinary MF linear scaling as in Eq. (|36|) . 
although this crossover disappears (e* 0) in the ther- 
modynamic limit. 

When e decreases further below e*, the system will 
reach the critical state where the particle density scales 
as pf ^ ^Jl/{gN) (see Eq. ^). The critical region in 
finite systems starts at e = Ec with 

Sc-y|-iV-V^withl/.= ^"(^^"^)/- , (39) 

where the finite size saturation starts to occur {2ec/g — 

The off-critical FSS behavior is illustrated in Fig.[T]^b). 
This scaling theory predicts that there exist two charac- 
teristic sizes TV* ~ e""* and Nc ~ e~'^ which separate 
three scaling regimes. In the regime I (e > £*) where 
N^, < N, the system behaves as in a SF network with 
infinite N and infinite fcmax, e.g. ps ~ e^. In the regime 
II (ec < £ < £*) where A'^c < A^ < A^*, it behaves as in 
a SF network with infinite A^ but with finite fcmax- The 
density scales as Ps — ^.ejg ^ sN^'^'^^'^'^/^ . Finally, it 
behaves as in a SF network with finite A^ and fcmax in the 
regime HI (e < £c) where A^ < Nc- The density scales as 

- A^-" with a = (l + (3-7)/w)/2 (see Eq. (PI) ). 

At the special case of the natural cutoff with lo — 7 — 1, 
the regime II disappears and there is a direct crossover 
from the regime I (no size effect) to the regime HI (critical 
size scafing). For 7 > 3 where 6 = i and g ~ 0(1), pf ~ 
e in both the regime I and II, so £* becomes meaningless. 
Here again we observe a direct crossover from the regime 
I to the regime HI. 

The FSS scaling behavior in the annealed SF network is 
sharply contrasted with that in the quenched SF network. 
While there are two characteristic sizes A^* and Nc that 
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FIG. 4: (Color online) (a) Scaling plots oi y — psN" vs. 
X = eiV^/" at 7 = 2.5 and to = 1.5, 2.0, 2.5, and 3.0. Network 
sizes are A'^ = 10"^, ■ ■ ■ , 10^. (b) Plots oi y/x against x. Each 
data set is shifted vertically by a constant factor to avoid an 
overlap. 



there is no flat region at = 1.5 (natural cutoff). As 
Lo increases, one can see clearly broadening of the flat 
region (regime II) which becomes larger with increasing 
N . This behavior is qualitatively consistent with the ex- 
pected FSS behavior. It is very difficult to observe the 
regime III scaling even at = 10^, similar to the diffi- 
culty encountered in the study of the critical dynamics 
in Sec. Hm 

Finally, the off-critical dynamic behavior can be easily 
derived from the rate equations, ([55)1 and ([M)) . in the 
thermodynamic limit, approaching the critical! ty from 
the active or the absorbing side: ps[t)—p^ ^ e^*/'^ where 
the relaxation time scales as r ~ £~''* with z^t = 1 in 
all cases. These results are consistent with our previous 
results through relaxation time relations of z = i^t/S' and 



SAMPLE-TO-SAMPLE FLUCTUATIONS 



depend explicitly on the degree cutoff in the former (at 
least for 2 < 7 < 3 and oj > 7— 1), it has been proposed in 
the latter through a droplet-excitation (hyperscaling) ar- 
gument Q that there exists a unique cutoff-independent 
characteristic size Nq ^ e^"" with 1/Dq — (7 — 2)/(7 — 1) 
for 2 < 7 < 3 and 1/Dq — 1/2 for 7 > 3. It is interesting 
to notice that the FSS theory in the annealed network 
coincides with that in the quenched network for 7 > 3 
and also at the special case with the natural cutoff with 
Lu — "f — I for 2 < 7 < 3. The origin of the discrepancy in 
the FSS theory between two different networks as well as 
the relevance/role of the quenched linking disorder have 
not been fully explored as yet, which awaits a further 
investigation. 

We have performed extensive simulations in the an- 
nealed SF networks to test the off-critical FSS the- 
ory. In Fig. m we present a scaling plot of ?; = PsN" 
against a scaling variable x = eiV^/" at 7 = 2.5 and 
Lu = 1.5,2.0,2.5, and 3.0. When w = 7 - 1 = 1.5 (nat- 
ural cutoff), the FSS theory predicts that the quantity 
y converges to a constant value for a: ^ 1 (regime III) 
and scales as ?/ ~ a;^/^'''"^) = a;^ for x ^ 1 (regime I). 
There should be no regime II. The numerical data in 
Fig. |4] seem to support this two-regime scaling behavior 
reasonably well. 

When w > 7 — 1, we expect three scaling regimes. Nu- 
merical data in the regime II and III will converge to 
a single curve, but those in the regime I should deviate 
from it because of the two different characteristic sizes. 
The numerical data in Fig. [U^a) show a clear evidence 
of the regime III (flat region), but a weak signature of 
the regimes II (linear-slope region, y ^ x) and I (no col- 
lapse). Although the signature is not prominent due to 
strong finite size effects, the existence of the three scaling 
regimes is evident. In Fig.dl^b), we present the numerical 
data in a different style by plotting y/x against x, so the 
regime II can be identified by a flat region. As expected. 



Suppose that one wants to generate a network of 
N nodes with a given degree distribution P{k) with a 
(forced or natural) cutoff. In general, there are two kinds 
of quenched disorder to be considered. First, one should 
sample a degree sequence {ki,--- ,A:jv} from P{k) and 
then choose a way of linking the nodes together to create 
a network. Disorder can be involved in both processes, 
which is named as sampling disorder and linking disorder, 
respectively. A quenched network involves both sampling 
and linking disorder, in general. 

An annealed network is free from the linking disorder. 
However, it may still have the sampling disorder. In the 
numerical studies in the preceding sections, we have sam- 
pled the degree sequence deterministically without any 
disorder. On the other hand, probabilistic sampling of 
the degree sequence leads to the sampling disorder. In 
this section, we investigate sample-to-sample fluctuations 
in annealed networks due to the sampling disorder. The 
quantity of our primary interest is g = (fc^)/(fc)^. 

When N values {ki, • • • , k^} are drawn probabilisti- 
cally in accordance with the distribution P{k) for fcmin < 
k < fcmax, a sampled distribution P{k) = '}2ii^i^k,ki/N 
may deviate from the target distribution P{k) due to the 
finiteness of N . The deviation is denoted by SP{k) = 
P{k) — P{k). Then, it is straightforward to show that 

[5P{k)] = (40) 

[5P{k)5P{k')] = -£Epn + ms,,„, (41) 

where [• • • ] denotes the sample (disorder) average. 

The nth moment of the degree of a sample is given by 

(fc") ^ J2 ^"^(^) = (1 + ^) ' (42) 

where we introduce shorthand notations as (n)o = 
^fe/fc"P(A:) and {n)s = 'Ek^'^^Pi^)- There is an 1/iV 
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factor in the correlator in Eq. (|4T|) . So, {n)s/ {n)o can be 
considered as a small expansion parameter for large N. 
Up to the second order, the quantity 5 of a given sample 
can be written as 



(2)o 



1 



(2)^ 
(2)o 



(l)o^ (1)2 (2)o(l>o 



The disorder-averaged correlators in Eqs. ((40|) and ((4T|) 
imply that [(jt.)^] = and that 

[{m)s{n)s] = ((to + n)o - (TO)o(n)o) . (43) 
This allows us to expand systematically [g] and (A^)^ = FIG. 5: (Color online) Probability distribution Pg{x) for x 




(.95 1 



.6 0.8 1 1.2 1.4 1.6 



[g^] — [g]'^ in powers of After some algebra, we obtain 



the following result up to the order 



[ 1 <2)o 



1 



o(2)o (3)o 
(l)o (l)o(2)o 



- 1 



(44) 



and 



iV(l)4 1(2)2 ^(i)o(2)o+ (1)2 ^ 

(45) 

Our interest lies in the SF network of N nodes having 
the degree distribution P{k) cx fc"''' in the interval fcmin < 
k < fcmax = N^/"^ with 7 > 2 and > 7 - 1. The 
term in Eq. (|44|) is always subleading, so the scahng 
behavior of [g] is determined by (2)o, which yields that 



- A^(3-7)/^ for 2 < 7 < 3 
[g] - <! ^logiV for 7 = 3 

- 0(1) for 7 > 3. 



(46) 



On the other hand, the term (4)o/ (2)§ in the parenthesis 
of Eq. (j45p makes a leading order contribution. Hence, we 
find that the relative variance Rg = (Ag)'^ /[g]'^ is given 
by 

^ for 2 < 7 < 3 

^ Ar2/»-i(iogiV)-2 for 7 = 3 
Ra = < ^ ^(5-7)/«'-i for 3 < 7 < 5 (47) 

- N-^logN for 7 = 5 

~- N^'^ for 7 > 5. 

In the theory of disordered systems, the relative vari- 
ance Rx of an observable X due to a quenched disorder 
is an indicator of the self-averaging property . When 
it vanishes in the thermodynamic limit N ^ 00, such 
a system is said to be self-averaging. The self-averaging 
property implies that an observable measured in a sam- 
ple with a typical disorder configuration takes the same 
value as the sample-averaged value in the N ^ 00 limit. 
A system with Rx ^ N^^ is said to be strongly self- 
averaging (SSA). This is the case when the central limit 



g/[g] in the annealed networks with A'^ = 10 
7 = 2.75 and to = 3.0 (b) 7 = 2.75 and to = 1.75. 



theorem works. When Rx ~ iV^*" with r < 1, such a sys- 
tem is said to be weakly self-averaging (WSA). A system 
with strong or relevant disorder lacks the self-averaging 
property near the criticality. In such a system, Rx con- 
verges to a finite value as increases. 

The result in Eq. (|47)) discloses the self-averaging prop- 
erty of the annealed SF network under the sampling dis- 
order. First of all, we find that the system with 7 > 5 
is SSA at all values of the degree cutoff exponent ui. 
For J < 5, Rg decays slower than at all values of 
w > 7 — 1. So the system is WSA. 

Interestingly, the systems lack the self-averaging prop- 
erty when 2 < 7 < 3 and w = 7 — 1 {Rg approaches a 
non-zero constant as N increases). Note that the cutoff 
exponent a; = 7 — 1 corresponds to the natural cutoff. 
Networks without explicit constraint on the degree also 
display the cutoff scaling /cmax ~ A^i/(t~i). In these net- 
works, not only the node-to-node degree fluctuation but 
also the sample-to-sample degree fluctuations are very 
strong. 

We present numerical data showing the (non-)self- 
averaging property in Fig. [51 Drawing N values of k from 
the distribution P{k) ^ A:"^ in the interval 2 < fc < A^", 
we calculated g = {k'^)/{k)^. This was repeated Ns — 
10^ times, from which one can construct the probabil- 
ity distribution function Pg{x) for x = g/[g]- Figure [5l^a) 
shows that the distribution becomes sharper and sharper 
as A^ increases. It indicates the self-averaging property 
at 7 = 2.75 and oj = 3.0. On the other hand. Fig. ^h) 
shows that the distribution converges to a limiting distri- 
bution. It indicates that the system is not self-averaging 
at 7 = 2.75 and w = 7 - 1 = 1.75. 

The strong disorder fluctuation raises an important 
question. In general, a real complex network is a dis- 
ordered media having quenched disorder, for example, 
the sampling disorder and the linking disorder as men- 
tioned before. Being coupled with dynamic degrees of 
freedom, the quenched structural disorder may give rise 
to disorder-relevant critical phenomena. This is a plausi- 
ble scenario, but has been ignored in most studies. It 
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seems a quite challenging problem to incorporate the 
quenched disorder into a systematic analysis. 

In the annealed network considered in this study, the 
dynamic degrees of freedom is completely decoupled with 
the sampling disorder (no linking disorder). Hence, the 
scaling theory developed here should be valid whether 
the sampling disorder is self-averaging or not. However, 
our result still warns that the sample average of any ob- 
servable involving g is practically meaningless due to its 
broad distribution, which occurs in the critical region 
(regime H and IH) in annealed networks with 2 < 7 < 3 
and the natural cutoff. 



VI. DISCUSSION AND SUMMARY 

We studied critical behavior of the CP in annealed 
scale-free networks. For the degree exponent 7 > 3, the 
standard single-parameter FSS is found with various dy- 
namic and static exponents which are independent of the 
cutoff exponent lo and also 7. For highly heterogeneous 
networks with 7 < 3, there exist two different character- 
istic time scales and their associated exponents depend 
not only on 7 but also on w. These results are con- 
trasted with those in quenched scale-free networks where 
a single-parameter FSS is found without any cutoff de- 
pendence even for 7 < 3 if the cutoff is not strong enough 
{uj < 7) d, . At the special case of = 7 — 1 (natural 
cutoff), these two different FSS coincide to each other. 

Annealed networks may include the sampling disor- 
der, which generates a strong sample-to-sample fluctu- 
ation in highly heterogeneous networks with the natu- 
ral cutoff. In quenched networks, the linking disorder 
is inherent, which generates the density-density correla- 
tion in neighboring nodes through coupling with fluctu- 
ating variables. This correlation leads to the shift of the 



transition point of the CP model fll, [T^ . In addition, 
the linking disorder generates another type of sample-to- 
sample fluctuations which cause spreading of the tran- 
sition points in finite systems. Hong, Ha, and Park [J 
showed that there exists a characteristic (droplet) size 
scale diverging as Ng ~ e""^' with l/Pq — (7— 2)/{j— 1) 
for 7 < 3. For N < Ng (or equivalently e < £q with 
Eq ^ 7V^(')'^2)/(7-i)-j^ ^Yie system feels the droplet length 
scale and the finite-size effect is dominant. As Sq > Sa 



given in Eq. (|39p . one may expect that the finite-size sat- 
uration induced by the droplet length scale comes in ear- 
lier (at £ = £q) \vl quenched networks than in annealed 
networks. Then, the cutoff-dependency of the satura- 
tion density may disappear. However, as £q < given 
in Eq. the cutoff-dependent density-decaying dy- 

namics comes in before saturation. The linking disorder 
fluctuation may be responsible for the disappearance of 
this dynamics in the quenched networks, but this is just 
a speculation as yet. A full understanding for the FSS 
behavior in quenched networks needs a further investiga- 
tion. 

During the final stage of preparing this manuscript, 
Bogufia, Castellano, and Pastor-Satorras posted a 
preprint on the cond-mat archive IJ|, the results of 



which partially overlap with those presented here. 



Acknowledgments 

This work was supported by KOSEF grant Accelera- 
tion Research (CNRC) (Grant No. R17-2007-073-0I001- 
0). This work was also supported by the Korea Research 
Foundation grant funded by MEST (Grant No. KRF- 
2006-003-C00I22). 



[1] S.N. Dorogovtsev, A.V. Goltsev, and J.F.F. Mendes, 

arXiv;0705:0010 (2007). 
[2] D.J. Watts and S.H. Strogatz, Nature 393, 440 (1998). 
[3] A.-L. Barabasi and R. Albert Science 286, 509 (1999). 
[4] H. Hong, M. Ha, and H. Park, Phys. Rev. Lett. 98, 

258701 (2007). 

[5] C. Castellano and R. Pastor-Satorras, Phys. Rev. Lett. 

100, 148701 (2008). 
[6] H. Hinrichsen, Adv.Phys. 49, 815 (2000). 
[7] K.-I. Goh, B. Kahng, and D. Kim, Phys. Rev. Lett. 87, 

278701 (2001). 

[8] M. Catanzaro, M. Bogufia, and R. Pastor-Satorras, Phys. 
Rev. E 71, 027103 (2005). 



[9] W.M. Hwang, S. Kwon, H. Park, and H. Park, Phys. 
Rev. E 57, 6438 (1998). 
[10] M.E. Fisher, J. Stat. Phys. 34, 667 (1984). 
[11] A. Aharony and A.B. Harris, Phys. Rev. Lett. 77, 3700 
(1996). 

[12] M. Ha, H. Hong, and H. Park, Phys. Rev. Lett. 98, 
029801 (2007). 

[13] C. Castellano and R. Pastor-Satorras, Phys. Rev. Lett. 

96, 038701 (2006). 
[14] M. Bogufia, C. Castellano, and R. Pastor-Satorras, Phys. 

Rev. E 79, 036110 (2009) re-print lajXiv:0810.3000p . 



